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We report measurements and calculations of the spin-subband depopulation, induced by a parallel 
magnetic field, of dilute GaAs two-dimensional (2D) hole systems. The results reveal that the shape 
of the confining potential dramatically affects the values of in-plane magnetic field at which the 
upper spin subband is depopulated. Most surprisingly, unlike 2D electron systems, the carrier- 
carrier interaction in 2D hole systems does not significantly enhance the spin susceptibility. We 
interpret our findings using a multipole expansion of the spin density matrix, and suggest that the 
suppression of the enhancement is related to the holes' band structure and effective spin j = 3/2. 

PACS numbers: 73.50.-h, 71.70.Ej, 73.43.Qt 



I. INTRODUCTION 



It was first pointed out by Janak that for a two- 
dimensional electron system (2DES) in a static mag- 
netic field the exchange interaction acts like an effec- 
tive magnetic field (in addition to the applied field) so 
that the Zeeman energy splitting is enhancedi Recently, 
the Zeeman splitting and spin susceptibility of inter- 
acting 2D carrier systems have been a subject of re- 
newed interest, 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 fueled by 
the promise of a paramagnetic to ferromagnetic ground 
state transition at very low densities; 19 ' 20 and the pos- 
sibility that the spin polarization is related to the ap- 
parent metal-insulator transition in dilute 2D systems. 21 
Experiments have mostly focused on determining the 
spin susceptibility from magneto-transport 
and magnetization^ measurements. The results gener- 
ally show that the spin susceptibility of 2DESs in dif- 
ferent materials, e.g., Si, 2 ' 3 ! 4 ! 5 ! 10 GaAs, 6 - 7 . and AlAs 8 - 9 . 
increases as the density is reduced, one report^ even sug- 
gesting a ferromagnetic instability at the lowest densities. 

Lately, the spin polarization of GaAs 2D hole sys- 
tems (2DHSs) has become the subject of intensive 
researchi 1 - 1 ^* 1 ^ 4 * 1 ^ 6 - 1 ^ because the holes have a larger 
effective mass (than electrons) so that they can be made 
effectively more dilute while maintaining high quality. 
Furthermore, the spin polarization of holes is important 
in the context of ferromagnetic semiconductors such as 
GaMnAs where it is known that the ferromagnetism is 
mediated by the itinerant valence band holesi^ 2 ^ We 
show here that the spin susceptibility of 2DHSs de- 
pends dramatically on the shape of the confining po- 
tential. Moreover, we find that, in contrast to their 2D 
electron counterparts, dilute 2DHSs exhibit no signifi- 
cant enhancement of the spin susceptibility as compared 
with calculations which neglect exchange-correlation^ 
We will argue that this surprising behavior is related to 
the holes' band structure and the fact they have effective 
spin j = 3/2 rather than j = 1/2 which is the case for 
electrons. 



TABLE I: Typical densities n (in 10 10 cm 2 ) and mobilities 
jU (in m 2 /Vs) of the samples used in this study. 
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II. SAMPLE PARAMETERS AND 
EXPERIMENTAL DETAILS 

Four samples from different wafers, including two 
GaAs/AlGaAs heterojunctions and two GaAs quantum 
wells (QWs) flanked by AlGaAs barriers, were investi- 
gated in this study (Table QJ. Depending on their sub- 
strate orientation and carrier type, our samples were ei- 
ther Be-doped (samples H, Ql) or Si-doped (Q2, A). 
All samples were fitted with metal front and back gates 
to control their density as well as the electric field per- 
pendicular to the 2D systems. We made measurements 
in 3 He or dilution refrigerators down to a temperature 
T = 0.03 K and in magnetic fields up to 25 T. 



III. EXPERIMENTAL RESULTS FOR (001) 2D 
HOLES 

In Fig. ^a) we show the longitudinal resistivity p xx 
versus in-plane magnetic field B» for samples H and Ql 
both measured at a density of n = 3.7 x 10 10 cm~ 2 . The 
data shows a positive magnetoresistance with a marked 
change in functional form above the magnetic field Bd 
that reflects the complete depopulation of the minority 
spin subband. 2-13 In Fig. ^a) Bd is marked by arrows. 
Remarkably, the field Bd depends greatly on the shape of 
the confining potential. Indeed, we have Bd — 10.6 T for 
sample H and Bd — 20.5 T for sample Ql, even though 
the data were taken at the same density. In Fig. ^c) we 
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FIG. 1: (a) Longitudinal resistivity p xx versus in-plane mag- 
netic field Bu measured at T — 0.3 K for 2D hole samples, H 
and Ql, at the same density n = 3.7 x 10 10 cm -2 . The depop- 
ulation fields Bd are marked by arrows, (b) Calculated density 
n_ in the minority spin subband of samples H and Ql as a 
function of Bu . (c) Measured (squares) and calculated (solid 
line) depopulation field Bd versus change A£ of the electric 
field in sample H for constant density n = 4.5 x 10 10 cm -2 . 
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FIG. 2: Measured (symbols) and calculated (lines) depopula- 
tion field Bd as a function of n for samples (a) H and (b) Q2. 
In (a), for the solid line and squares (dashed line and circles) 
n was varied via a back (front) gate. In (b) the different sym- 
bols refer to B || [TlO] and B || [332] as indicated. The upper 
horizontal axes show the calculated density paramater r s for 
the corresponding n. In (a) we used r s (n) for the front gate. 

show Bd in sample H when the electric field £ across the 
junction is varied by means of front and back gates such 
that n is kept constant at 4.5 x 10 10 cm~ 2 . The field Bd 
increases significantly with increasing £. In Fig. [21 a) we 
show the measured Bd versus n for sample H. The values 
of Bd depend rather sensitively on whether n is changed 
by means of a front or back gate. 

IV. CALCULATIONS AND DISCUSSIONS 

In order to explain the experimental results of Figs.^ 
and [3] we have performed parameter-free calculations 
in the multiband envelope-function and self-consistent 
Hartree approximations for the quasi-2D systemic Fig- 
ure n^b) shows the calculated density n_ in the mi- 



nority spin subband as a function of Bu. The lines in 
Figs. Ulc) and [3 show the calculated Bd for the corre- 
sponding experiments^ The calculations reproduce the 
different behavior of samples H and Ql in satisfactory 
agreement with experiment. 



A. Confinement potential dependence of Bd 

As we discuss in the next section, the close agree- 
ment between the experimental and calculated Bd in 
Figs, n and [^a) is very surprising because the calcu- 
lations do not take exchange-correlation effects into ac- 
count. Such effects are indeed dominant for 2D electron 
systema 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 18 ! 27 ! 28 ! 29 that are as dilute as the 
2DHSs of Figs.n an d|2I a )- Before elaborating on this as- 
pect of our results, however, we first discuss the remark- 
ably strong dependence of Bd on the confining potential. 
For an ideal, strictly 2D system with effective mass m* 
and effective g factor g* we have Bd cx l/(m*g*) in- 
dependent of the shape of the confining potential. To 
understand the surprising results in Figs.n an( i0we will 
first concentrate on the Zeeman splitting which gives rise 
to the dominant contribution of the confinement depen- 
dence of Bd in 2DHSs. Then we discuss the effect of Bu 
on the orbital motion. 

Unlike electrons in the conduction band that have spin 
1/2, holes in the uppermost valence band are character- 
ized by an effective spin 3/2 (Ref. |25T> . Subband quanti- 
zation in 2DHSs yields a quantization of angular mo- 
mentum with z component m = ±3/2 for the heavy 
holes (HHs) and m = ±1/2 for the light holes (LHs). 
In our samples, only the lowest HH subband is occu- 
pied. The quantization axis of angular momentum that 
is enforced by HH-LH splitting points perpendicular to 
the 2D plane. The Zeeman energy splitting due to Bu 
thus competes with the HH-LH splitting and it is well- 
known that the By-linear Zeeman splitting of HH states 
is suppressedi 30 i 31 [The simple model of Ref. [3(J yields 
Bd — 250 T for the systems in Fig. Ufa).] I n the follow- 
ing we will discuss why the depopulation fields Bd ob- 
served in real 2DHSs are much smaller than what these 
arguments suggest. 

The dispersion of HH states is known to be highly non- 
parabolic as a consequence of HH-LH coupling. 25 There- 
fore, the suppression of Zeeman splitting linear in Bu is 
merely the lowest-order effect in a Taylor expansion of 
the spin-split dispersion ^(feii, Bu) of HH states as a 
function of the (canonical) wave vector fey , Bu , and spin 
index a. Mixed higher-order terms proportional to Bu 
and feii give rise to an average Zeeman splitting of the 
occupied hole states which is approximately linear in B\\ . 
Thus we find that Bd is generally much smaller than the 
value one would expect if the Bu -linear Zeeman splitting 
were suppressed. This is also consistent with previous 
experimental data for 2DHSs that were interpreted ig- 
noring completely the suppression of Bu -linear Zeeman 
splitting in HH systems. 11 ' 12 ^ 3 ! 14 ! 15 
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Now we can understand why the Zeeman energy split- 
ting in 2DHSs depends sensitively on the shape of the 
confining potential. The mixed higher-order terms that 
are responsible for the Zeeman energy splitting Ez(B^) 
of HH systems compete with the HH-LH splitting. The 
latter depends sensitively on the shape of the confining 
potential so that we have here a tool to tune Ez(B^) of 
2DHSs. In narrow quasi-2D HH systems we have a large 
HH-LH splitting so that the Zeeman energy splitting is 
reduced, giving rise to a large B d . We get a large E Z {B\\) 
(a small B d ) in wide systems. We can define E z {B d ) 
as the energy difference between the Fermi energy and 
the subband edge at B d . In the wide heterojunction of 
Fig. [lfb), the calculated E z (B d ) is 0.44 meV, signifi- 
cantly larger than Ez(B d ) = 0.26 meV in the narrower 
QW, despite the smaller value of B d in the heterojunc- 
tion. Similarly, the increase of B d with increasing AS in 
Fig. ^c) reflects the change of the HH-LH splitting in 
the system. 

Next we discuss the effect of B\\ on the orbital mo- 
tion. In general^ the mass of the particles in quasi-2D 
systems increases as a function of B\\ which reflects the 
fact that, ultimately, for large B\\ resulting in a mag- 
netic length comparable to the width of the quasi-2D 
system the particle states become dispersionless Landau 
levels. Obviously, this effect depends on the thickness 
of the quasi-2D system and it has been shown that B d 
in wide quasi-2D electron systems is much smaller than 
B d in narrow 2DESsi£ We will argue next that the mass 
enhancement does not explain, however, the results in 
Figs. [Hand Ufa). 

Our numerical calculations show, in agreement with 
the 2DESs' results; 7 - that the mass enhancement at small 
Bii is smaller in the QW than in the heterojunction. 
However, m* in 2DHSs increases highly nonlinearly as 
a function of Bit which is particularly important for the 
QW with the larger B d . Thus we find that at B d the 
mass enhancement in the narrower 2DHS of the QW is 
larger than in the wide 2DHS of the heterojunction. We 
note that at B d the mean kinetic energy equals approxi- 
mately half the Zeeman energy splitting Ez{B d ) so that 
for Fig.^b) the mass enhancement can be inferred from 
the Ez values quoted above [see also Eq. below]. 

The anomalous enhancement of m* at B d with decreas- 
ing width of the quasi-2D HH system depends sensitively 
on the system parameters such as the density and the 
shape of the confining potential. For the paramteres in 
Fig. W[p), m* at B d is approximately independent of AS 
(despite the significant change of B d ), i.e., the increase 
of B d with AS is essentially only due to the decrease of 
the Zeeman splitting discussed above. For about twice 
the largest field AS we could reach experimentally one 
enters the regime when m* at B d starts to increase with 
AS. 
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FIG. 3: Ratio B%/B^ p of the depopulation field B° d calcu- 
lated neglecting exchange-correlation to the measured field 
B7 P for a 2DES (squares) and a 2DHS (circles) in GaAs (001) 
heterojunctions plotted versus the density parameter r a . In 
both cases, n was varied via a front gate. For electrons, the 
solid line shows the ratio B®/Bj c , where B^ c was calculated 
taking into account exchange-correlation^*— ■ 



B. Lack of spin susceptibility enhancement for 
(001) 2D holes 

A most remarkable aspect of the results in Figs. 
and EJa) is the reasonable quantitative agreement be- 
tween the experimental data and the calculations. 
This is particularly puzzling because many-particle ef- 
fects beyond the Hartree approximation (i.e., exchange- 
correlation effects) were not taken into account. This 
is in sharp contrast to the case of dilute 2DESs 
for which it is known that exchange-correlation sig- 
nificantly increases the spin susceptibility when n is 
reduced, 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 18 ! 27 ! 28 ! 29 To quantify this point 
we show in Fig. |3| the ratio B^/B^ of the depopula- 
tion field B®, calculated neglecting exchange-correlation, 
to the experimentally measured field B d xp for a 2DES 
(squares) 33 and the 2DHS (circles) in sample H. The ratio 
thus reflects the enhancement of the spin susceptibility at 
B d due to exchange-correlation. Our results are plotted 
as a function of the dimensionless density parameter r s 
defined as the average interparticle spacing measured in 
units of the effective Bohr radius tig, r s = l/(a^y / nn). 

For the 2DES in Fig. the ratio B° d /B° d xp is be- 
tween 2 and 4 and, as expected, it increases with 
r s . We also remark that for electrons the experimen- 
tally observed reduction of B d is in reasonable quanti- 
tative agreement with numerical calculations that take 
exchange-correlation into account. To illustrate this 
point, the solid line in Fig. [3] shows the ratio B d /B d c , 
where B d c was calculated in the framework of spin- 
density-functional theory using a parameterization of 
the polarization-dependent exchange-correlation poten- 
tial that was recently obtained by means of quantum 
Monte Carlo calculationsi 7 * 2 ^ For the 2DHS, on the other 
hand, the expected enhancement of the spin susceptibil- 
ity and B d /B^ p ratio is conspicuously absent in Fig.|3] 
Note that, because of their larger effective mass com- 
pared to GaAs electrons (m* ~ 0.25 compared to m* = 
0.067; here m* is given in units of the free-electron mass), 
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2DHSs have significantly larger r s and are thus effectively 
much more dilute. Nonetheless, the ratio B^/B^ xp re- 
mains close to unity up to the largest values of r s where 
a greater than ten-fold enhancement is expected. 

Before discussing possible reasons for this anomalous 
behavior of 2DHSs, we make remarks regarding the effec- 
tive mass to* which enters a|j and thus determines r s . For 
holes, to* is not uniquely defined. As discussed above, 
the HH dispersions are typically nonparabolic, meaning 
that to* depends on energy and therefore on n and the 
confinement potential. Moreover, the HH systems have 
a lar ge R ashba and Dresselhaus spin splitting at B = 
fRef. |25T) . leading to two energy versus wavevector (fen) 
dispersions with different curvatures and effective masses 
ml and m*_. Commonly, values of m* between about 0.2 
and 0.4 are used for holes in GaAs»iA*i^i£i2&i2i Here we 
adopt a simple definition for an average effective mass 
(to*): 



(m*) = h 2 TTn/(2(E k }) 



(1) 



where (Ek) is the mean kinetic energy per particle. Fig- 
ure Ufa) shows the calculated density parameter r s in 
sample H, when n is changed by means of a front or 
back gate. Note that for a single, parabolic dispersion 
with an effective mass to*, the mass (m*) as defined in 
Eq. m properly reduces to to* and is independent of 
n. For the 2DHS, on the other hand, (m*) in general 
depends sensitively not only on n but also on the sys- 
tem's parameters such as the thickness of the 2DHS and 
on the applied electric and magnetic fields, as discussed 
above. If we take into account Rashba and Dresselhaus 
spin splitting, 25 then we get, similar to Eq. effective 
masses (toj.) for each spin subband. To illustrate this 
effect, we show (toj_) for sample H in Fig.^Jb). We em- 
phasize that the main conclusion of our work, namely 
the lack of enhancement of the spin susceptibility with 
increasing diluteness, is not affected by the specific val- 
ues of to* used to define r s : it is clear in Fig. [3] that if r s 
were changed by a factor of 2 or 3, there would still exist 
a large discrepancy between the experimental hole data 
and the expected enhancement. 

Why do dilute 2DHSs not show a significant enhance- 
ment of the spin susceptibility? Using a recently devel- 
oped multipole expansion of the spin density matrb«i& we 
argue in the remainder of this paper that the j — 3/2 hole 
spin is the likely culprit. For 2DESs with spin 1/2 it is 
well-known that the mean Coulomb energy (E c ) per par- 
ticle can be completely characterized using n and (the 
magnitude of) the spin polarization £ as independent 
parameters^ This is because the 2x2 spin density ma- 
trix of spin 1/2 systems can be decomposed into four 
independent terms: n (a monopole) and the three com- 
ponents of the spin polarization vector £ (a dipole) ^ In 
the Hartree-Fock (HF) approximation, the direct part of 
the Coulomb energy (E^ F ) cancels the potential of the 
positive background so that only the exchange term (E x ) 
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FIG. 4: Calculated density parameter r a , average effective 
mass (m±), normalized dipole and octupole moments— versus 
n for samples H and Q2. Line styles have the same meaning 
as in Fig. [5] 
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(E?*) = (E X ) = 



(i + c) 3/2 + (i - C) 3/2 



The Coulomb energy (E^ F ) is thus proportional to s/n. 
Higher order terms in a series expansion for (E c ) (n, C) of 
2DESs were calculated in Ref. \27\ (E c ) was calculated 
numerically in, e.g., Refs.l28l andl29l 

The 4x4 spin density matrix of j = 3/2 2DHSs, on 
the other hand, can be decomposed into four multipoles, 
where the monopole is the density n, the dipole corre- 
sponds to the spin polarization at B > 0, the quadrupole 
reflects the HH-LH splitting, and the octupole is a unique 
feature of j = 3/2 hole systems at B > (Ref. Il6|) . For 
sample H, the normalized dipole and octupole at Bd are 
shown in Figs. He) and (d) (Ref.©. Unlike 2DESs, the 
dipole at Bd is much smaller than unity, i.e., despite the 
fact that only one spin subband is occupied at Bd, the 
system is only weakly spin polarized. This resulliii is an 
immediate consequence of the suppression of Bu -linear 
Zeeman splitting22i2i discussed above. The octupole can 
be interpreted as a new "spin degree of freedom" of spin 
3/2 hole systems at B > which does not exist for the 
more familiar case of spin 1/2 electron systems. When 
the spin polarization is suppressed for an in-plane mag- 
netic field, the 2D HH systems aquire instead a large oc- 
tupole moment^ as visible in Fig.^Jd). The quadrupole 
is always close to unity because the HH-LH mixing is 
small in the systems considered here. Therefore, the 
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quadrupole is not shown in Fig. By definition, these 
four multipolcs provide a set of independent parameters 
that can be used to parameterize the Coulomb energy 
(E c ) of spin 3/2 systems, similar to (E c )(n, £) in spin 1/2 
systems. - 9 However, the series expansion is presently not 
known and its calculation represents a formidable task. 
Our study indicates that the series expansion of (E c ) of 
spin 3/2 2DHSs is qualitatively different from (E c )(n,C) 
of spin 1/2 2DESs. 

The HF exchange energy (E x ) of 2D HH systems at 
B = is the same as (E x ) of spin 1/2 2DESs because 
Eq. J2)l requires only that the eigenstates of the two spin 
subbands for the same fen are orthogonal^ For a HH 
system, the main effect of a perpendicular magnetic field 
B^ is a spin polarization (a dipole) , whereas an in-plane 
field Bu usually gives rise to an octupole moment>iSi2£ 
The spin density matrices of 2D HH systems at B± > 
and Bii > are thus qualitatively different. However, the 
HF exchange energy does not distinguish between these 
cases and always leads to the same enhancement of the 
exchange energy as in 2DESs. We note that different re- 
sults can be obtained for (E x ) when HH-LH mixing is 
significant 42, Also, different results are obtained for (E c ) 
in higher-order perturbation theory when the more com- 
plicated energy dispersion must be taken into account. 
These are the reasons why the well-established results 
for exchange-correlation in dilute spin 1/2 2DESs cannot 
easily be transferred to spin 3/2 2DHSs. 

V. RESULTS FOR (113) 2D HOLES 

We extend our investigation by comparing the results 
for sample H with the data for Q2, a QW grown on a 
(113)A GaAs substrate. Figure Hfb) shows Bd versus n 
for Q2. The field Bd strongly depends on whether Bu is 
applied in the in-plane crystallographic directions [110] 
or [332] (Ref-El). The right column of Fig. H shows r s , 
(m5_), the dipole and the octupole moments calculated 
for Q2. For this sample, when By is applied parallel to 
[332], the measured Bd is well below the calculated value. 
It is remarkable that for this particular geometry, the 
octupole remains small, but the 2DHS develops a large 
dipole moment [Figs.EJc') and (d')], similar to 2DESs in 
a Bu CRef. li^l . This observation suggests that the spin 
susceptibility is enhanced by many-particle effects only 
when the magnetic field gives rise to a spin polarization. 
On the other hand, Figs. El and 01 suggest that there is no 



significant enhancement in j = 3/2 2DHSs with a large 
octupole but a small dipole moment. 



VI. SUMMARY AND CONCLUSIONS 

We have shown that the spin susceptebility of dilute 
GaAs 2DHSs in an in-plane magnetic field Bu depends 
sensitively on the shape of the confining potential. Most 
remarkably, the spin susceptibility is not significantly en- 
hanced as compared with calculations which neglect the 
carrier-carrier interaction. This is in sharp contrast to 
dilute electron systems for which it is known that many- 
body effects greatly enhance the spin susceptibilty. Using 
a multipolc expansion of the spin density matrix we have 
argued that the suppression of the enhancement is related 
to the holes' band structure and effective spin j = 3/2. 

Our findings have important implications for the quan- 
tum phase diagram of dilute 2DHSs. In dilute electron 
systems, the exchange-correlation enhancement of the 
spin susceptibility can be considered a precursor for the 
ferromagnetic liquid which is expected to be the ground 
state of ultra- low density 2DESs with r s > 26 (Ref.|H). 
The extra multipolcs of 2DHSs provide new possibilities 
for the ground state of hole systems to respond to ex- 
ternal perturbations such as a magnetic field thus lead- 
ing to a richer phase diagram than in spin 1/2 elec- 
tron systems^ However, our results suggest that a fer- 
romagnetic phase (i.e., a fully spin-polarized phase with 
a maximum dipole moment) is often not favored in dilute 
2DHSs. This could also have important implications for 
ferromagnetic semiconductors such as GaMnAs where it 
is known that the ferromagnetism is mediated by the itin- 
erant spin 3/2 holes in the valence bandi^SiSi In itinerant 
ferromagnets it is the polarization-dependent competi- 
tion between the Coulomb energy and the kinetic energy 
of the interacting carriers which controls the ferromag- 
netic transition. 
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